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Abstract 

A kinetic model of gravimagnetic shock waves (GMSW) in a locally 
<^> anisotropic plasma is investigated. The equations of a drift approximation 

are written, and the moments of the distribution function are calculated. 
(— i ■ Solutions of the drift equations for a highly anisotropic ultrarelativistic 

plasma are found. It is shown that in this case the GMSW essentially 
affect the angular characteristics and the intensity of the magneto - brem- 
■ sstrahlung of the magnetoactive plasma. 

^ 1 Introduction 

£~ ' i In Ref. pp, from the requirement that the dynamical velocity of plasma be 

equal to that of the electromagnetic field due to the Einstein equations and the 
first group of the Maxwell equations, the equations of the relativistic magnetic 
hydrodynamics (RMHD) of magnetoactive plasma in a gravitational field were 
derived. In the same work (see also [2]) a remarkable class of exact solutions 
J ' of the RMHD equations was obtained; it describes the motion of a magnetoac- 

tive locally isotropic plasma in the field of a plane gravitational wave (PGW) 
and is called gravimagnetic shock waves (GMSW). In Ref. [3J it was shown 
that in pulsar magnetospheres the GMSW are a highly effective detector of the 
gravitaional radiation of neutron stars. An observed consequence of the energy 
transformation from the gravitaional wave to the GMSW energy are the so- 
called giant pulses, sporadically arising in the radiation of a number of pulsars. 
The estimates made in [3J and [3] allow one to identify the giant pulses in the 
pulsar NP 0532 radiation with the gravitational radiation of this pulsar in the 
basic quadrupole mode of a neutron star. 

A fundamental importance of GMSW for the gravitational theory leads to 
the necessity of a more detailed and comprehensive investigaton of this phenom- 
ena. As was shown in pQ and [3] , a GMSW is realized in an almost collisionless 
and nonequilibrium plasma situated in an abnormally strong magnetic field. 
Under these conditions, as a consequence of strong magneto - bremsstrahlung, 
the isotropy of the local distribution of plasma electrons is significantly vio- 
lated, as was assumed in obtaining the solution in [T]. In [5], on the basis of the 
general RMHD equations, a hydrodynamic model of GMSW in an anisotropic 
plasma was built. Since in the hydrodynamical approach the number of equa- 
tions obtained is smaller than that of unknown functions, we had to postulate 
a relation between the longitudinal and transverse components of the plasma 
pressure. In [5] the simpliest variant of such a (linear) relation was studied. 
The study revealed a high dependence of the GMSW process on the degree of 



X 

5-H 



1 



plasma anisotropy, leading to the necessity of building a dynamical model of 
anisotropic magnetoactive plasma motion in the gravitaional radiation field. It 
is the problem dealt with in the present paper. Throughout the paper a set of 
units is used where (c = G = h = 1). 



2 Drift solution to the kinetic equation for an 
anisotropic plasma 

2.1 Field quantities in the plane gravitaional wave metric 

Let us study a collisionless plasma in the PGW metric [5]: 

ds 2 = 2dudv-L 2 [e 2fi {dx 2 ) 2 +e- 2fl (dx 3 ) 2 } = 2dudv-A(dx 2 ) 2 -B(dx 3 ) 2 , (1) 

where /3(u) is an arbitrary function (the PGW amplitude), the function L(u) 
(the background factor of the PGW) obeys a second-order ordinary differential 
equation; u — (t — a; 1 ) is the retarded time, v = ^= (i + x 1 ) is the advanced 
time. The absolute future corresponds to the range T + : {u > 0;v > 0}, the 
absolute past to T~ : {u < 0; v < 0}. The metric (QJ admits a group of motions 
£/5, with the corresponding three linearly independent Killing vectors at a point: 

e=si; e=si; e=si. (2) 

(1) (2) (3) 

Let there be no GW for u ^ 0, i.e., - 

P(u)\u^q = 0; L(u)\ u ^ Q = l, (3) 
and the homogeneous magnetic field be directed in the plane {x 1 ^ 2 }: 

#l|u<0 = H cos fi; H 2 \ u ^o = H Q sin Q,; 

^3|«^o = 0; ^i|u<0 = 0' (4) 

where is an angle between the Ox 1 axis (PGW propagation direction) and 
the magnetic field direction H. The conditions (@1 correspond to the vector 
potential 

A v = A u = A 2 = 0; A 3 = H^x 1 sin 17 - x 2 cosfi); (u < 0). (5) 

The electromagnetic field in the PGW metric ([T]) for an initially homogeneous 
plasma is described by the vector potential pQ: 

A 2 = A v = A u = ; 

^3 = — Hqx 2 cos Q H — -=Hq[v — ip{u)] sinO, (6) 

where tp(u) is an arbitrary diffcrentiable function satisfying the initial condition: 

VV<0 = (7) 
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In this case the only nonzero component of the Maxwell tensor depending on ip 
is 

Fua = —j=Hoil>' sin SI. (8) 
Other nontrivial components of the Maxwell tensor are: 

F 23 - -H cos SI; F v3 = -^=H sin St . (9) 

2.2 Collisionless kinetic equation 

The collisionless kinetic equation for the 8-dimensional distribution function for 
charged particles of a kind a, F a (x l ,Vj) has the form [7]: 

where 

Vi=p i + e a A i (11) 
is the generalized momentum of a particle and 

n a {x\Vi) = X -g l i (Pi - e a Ai) (V 3 - e a Aj) (12) 

is the Hamiltonian function of a charged particle. 

The process of obtaining equations in the drift approximation on the basis of 
the collisionless kinetic equations (fTUl was described by the author in [3]. Here 
we only somewhat transform the process for the case of an initially anisotropic 
distribution. Note, besides, that in the case SI ^ ir/2 the results [8] are erroneous 
because they use V% in the role of an integral of motion, which is the case only 
if SI = ir/2. That was mentioned in the work cited above pQ. 

Let us study the case when the GW propargates perpendicularly to the 
magnetic field, SI — ir/2, and let us look for solutions of the kinetic equation 
independent of the variables v,x 2 and x 3 . As in [8], let us introduce the unit 
timelike vector v l : 

v 2 =v 3 = 0; v u = ^; v v = (13) 

and transform the kinetic equation into the frame of reference (FR) moving with 
the velocity v l : 

V v = v u {Ti + Vx ) V u = Vv {Va - Vi ) . (14) 
The Jacobian of this transformation is equal to unity: 

w^) =1 - (15) 
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As a result, we get the equation: 

Cry , -77 \ dFa f f<n -n ^ (n dF * , 9F a ^\ ei? -n 



+ \ u-'yn + (B-'yvi] Vu + = o, (i6) 

where: 

2 e 2 - 9 
H 2 = --d u A 3 d v A 3 = H^'— (17) 

is an invariant of the electromagnetic field (the magnetic field intensity squared 
in the FR moving with the velocity v z ). 

On the PGW front the solution (|T5|) must satisfy the initial condition cor- 
responding to a homogeneous anisotropic currentless plasma: 

F a {x\Vi)\ u = = f a (pi,pl)5(U a - l -ml) , (18) 

where f a is an arbitrary function of its arguments and the following notations 
are introduced: 

P± =Pi +pI ; P\\ = -p-2 ■ ( 19 ) 
The kinetic equation (fTo) has three exact integrals: 

U a {x l ,Vi) = ^m 2 , = Const; 

Vi = Const ; V 3 = Const . (20) 
For a complete solution of the problem we need one more independent integral. 

2.3 Drift approximation 

We shall solve Eq. (|16p in the drift approximation, when the Larmor frequency 
for each kind of charged particles - 

UJ a = (21) 

m a 

is much greater than the characteristic frequency oj of the gravitational wave: 
A = — « 1 . (22) 

In the zeroth order with respect to the parameter A (Tl6]) takes the form: 

H wr^ (23) 



4 



i.e. in the drift approximation F a is independent obviously of V\. Thus in the 
drift approximation, apart from the above exact integrals, the kinetic equation 
has also drift (approximate) integrals [5]: 



The solution of the kinetic equation corresponding to the drift approximation, 
which, in the absence of GW, is transformed to an anisotropic and currentless 
one, can be written in the form 



where f a is an arbitrary function of its arguments. In particular, the following 
fa can be chosen: 



FR (fTi|) moving with the velocity vi the drift solution locally coincides with the 
unperturbed distribution (IT51) . 

Substituting the obtained distribution function of the zero drift approxima- 
tion (l2~5j) into the kinetic equation (ITB1) . it is easy to derive a correction of the 
first drift approximation, 5F a ~ A , for the distribution function determined 
up to an additive component being an arbitrary function of the above drift inte- 
grals (see [8]). However, exact consequences of the collisionless kinetic equation 
IT6|) are conservation laws for the number of each kind of particles and the total 
energy-momentum tensor (EMT) of the plasma and the electromagnetic field 
[5]. These laws impose certain restrictions on the above additive component 
and lead to differential equations for the functions H\\(u) and fj.±(u). 

Due to the symmetry properties, it turns out that the first order correction 
to the distribution function contributes only to the component 713 of the particle 
number current density vector and the components T" 3 and T„ 3 of the particle 
EMT. However, the component n3(u, v) does not affect the continuity equation, 
only the above two components of the particle EMT appear in the transport 
equation of the total EMT for the component i = 3 they do not appear in other 
transport equations. On the other hand, the drift current determined by the 
first correction to the distribution function can be obtained as a consequence of 
the energy- momentum conservation laws and the Maxwell equations, without 
addressing to a solution of the kinetic equation (see [1]). As a result, it turns 
out that the set of the transport equations and the Maxwell equations split 
into two subsets, one of which, determined by the zero drift approximation, is 
self-consistent and closed and entirely determines the functions ip(u), /J.\\(u) and 
/j,±(u). The correction of the first drift aproximation to the distribution function 
does not affect the magnetic hydrodynamics equation. 



V4 « Const ; u « Const . 



(24) 




(25) 
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3 Derivation of the magnetic hydrodynamics equa- 
tions 

3.1 Algebraic structure of the distribution function 
moments 

Returning to the original FR by means of (ThO) . let us introduce the following 
scalars: 

P\\ = (p, n); p 2 = (uV - g lj )PiPj ; 

P± = P 2 - p\ = (v l v j - g lJ - n l n :i )p. l p J , (27) 

where n 1 is the unit spacelike vector in the direction of the magnetic field in- 

* 

tensity vector Hi = v 3 Fji- 

= H 2 = -(H,H). (28) 

then the distribution function in the zero drift approximation (|26[) takes the 
form 

/a=/ (^ 2 pf + M: 2 pi)- (29) 
It is not difficult to show that the moments of this ditribution are: 

n\{x) = J f a (x,p)p l dP; 

P(x) 

T i J{x)= J .f a {x,p) P ydP, (30) 

P(x) 

where 

f=g(2S + l)dp 1 dp 2 dp 3 



dP - 

(S - is the particle spin), and have the following algebraic structure: 

K(u)=nv i ; (31) 

T l J{u) = (e+ P±yv j - P ± g ij + (P\\ - P ± )n l n j , (32) 

where n, e, Py and P± are some scalar functions of the variable u; the indices 
of the kind of particles, a, are dropped in these scalars for the sake of simplicity 
of notations. The EMT track of particles is equal to: 

r o («)=e-2Pj.-P|.. (33) 
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From (pTTj) it follows that the velocity vector vi introduced in ([TB")) coincides with 
the plasma kinematic velocity vector. It is not difficult to check the fulfilment 
of the condition [T] : 

(«,n) = 0, (34) 

and consequently the vector v % is an eigenvector of the particle EMT, i.e. it is 
the dynamic plasma velocity vector according to Synge |10j . thus the frame of 
reference introduced by the relations (|14j) is comoving the plasma. 

Note that in the first drift approximation the equation for kinematic and 
dynamic velocities of particles is violated and therefore a drift current arises. 

3.2 Magnetic hydrodynamics equation for an anisotropic 
plasma in the PGW field 

In Sec. 3 of Ref. p] the RMHD equations for a plasma in an arbitrary gravi- 

v 

tational field and for an arbitrary structure of the plasma EMT, Tij, with the 
eigenvector v % were obtained: 

Tij v 1 = ev t , (35) 

where the invariant e > is the plasma energy density in the comeoving FR. 
In this case the following conditions were imposed on the invariants of the 
electromagnetic field: 

/••'/•,, 0: (36) 
— > . (37) 
In [1] it was shown that, under these conditions, from the conservation law 

/;,' = (38) 
for the whole EMT of the plasma and the electromagnetic field 

p f 

ij ~Tij + Tij 

and the first group of the Maxwell equations 

" (39) 

it follows: 

1. The conditions of embedding of the magnetic field in the plasma are: 

FijV j = (40) 

(in this case the velocity vector vi also automatically becomes an eigen- 
vector of the electromagnetic field EMT) ; 
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2. The second group of the Maxwell equations is 

F% = -AnJ^ (41) 

with the drift current 

p 

J dr = F . F j m > ( 42 ) 

3. The differential relations are 

v ir r\ k =0, (43) 



iTT\ fc =0. (44) 

It is not difficult to ascertain that the Maxwell tensor determined by the 
vector potential ([6]) automatically satisfies the conditions (|37| and the velocity 
vector [T3"|) - the embedding conditions (j40|) . Therefore in this case, due to the 
EMT conservation (|39|) . which is an exact consequence of the kinetic equations, 
the relations (|41|) — (|44|) must hold. The isotropic Killing vector ([2]) gives an 
exact integral of Eqs. ([55]): (|3"9"1) : 

L 2 T^ = Const . (45) 

thus, taking into account Eqs. (fT3"|) and (|32p and the initial conditions ([3]) and 
(JTJ), we get from (T4"5|) : 

L 2 ( £ + P ± ) = V/(? + P±)A( U ). (46) 
where we have introduced the so-called GM^M^ governing function (see pQ): 

A(u) = l-a 2 (e 213 -1) (47) 
and the GMSW dimensionless parameter 

* 2 = f\ ■ (48) 
4ir(e + P ± ) 

It can be further shown that the relation (|44|) is transformed into an identity, 
and Eq. (g5J) gives: 



£ ' + (lnff)'(Pn-P ± )-^ln^ (e + Pj|) = 0. (49) 

The EMT conservation equation and the Maxwell equations do not give other 
independent relations. 
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The missing equations are obtained from the particle number conservation 
laws which are also an exact consequence of the kinetic equations: 



Thus we get the set of equations (|46|) . (|49|) and (ISTJl) for determining the 
three unknown scalar functions ip, /iy and It is only necessary to obtain in 
this way explicitly the scalars n, e, Py and P± out of (1 + 2n) scalars, where n 
is the particle kind number. 

Note that there are at least two kinds of charged particles (in the case of 
interest these are protons and electrons) in the plasma. Thus, there are two 
kinds of scalars /jy and fj,±, and one particle number conservation law for a given 
kind of particles connects each couple of them ([ST)]) . The summed components of 
pressure and energy density appear in equations (1461) and (1491) . As a consequence 
of the initial electroneutrality of the plasma, (JSHJ) , we get a relation between the 
local concentrations of neutrons and protons [8] : 

n e (u) = n p (u) . (51) 

3.3 Calculating the moments of the distribution function 

For calculating the above scalar functions it is necessary to find the moments 
of the distribution function (|30|) . An easiest way to do it is to use Eqs. (|3 1 1) 
and (|32|) in the comoving FR according to (Q~4|), using the property (fl~5]) of this 
transformation and using the spherical coordinates in the momentum space: 

Pi = [i\jp cos 9 cos </> ; 

Pa = (J-±P cos 9 sin 4> ; (52) 
Vi — /iyp sin 

In so doing we obtain: 



n\ = L- 2 (L 2 nv v )' = => 




(50) 



(53) 



where 



oo 




(54) 



o 



the integration variable is: 



2 -2 2 i -2 2 

P = V\\ P\\ + M_L P_L • 

thus we obtain from ([53]), ([50]) and (j5i~]) ): 



(55) 




(56) 
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Further, setting for definiteness 

M-L < M|| ) 

we get the expressions for the components of the plasma pressure and energy 
density: 



(25 + 1)m1m| f o 



f {v )v dv x 



^Jm 2 + fi 2 p 



2Z2 m + VIP 



In 



77l 2 + /in|5 2 + P.. //in - 



\/ 771 2 + fi\p 2 



;(57) 



P± = 



00 



r 1 

/ (p 2 )p 2 dp x — Jm 2 + fi 2 p 2 



m 2 + p 2 (2fi 2 - fi\) 



In 



mf + fif,p z + py hir, - fij_ 



Pa/Mm - H± 



a/to 2 + fi^p 2 



(58) 



e = 



25+ 1 
4tt 2 



«_M|| / f(p 2 )p 2 dpx 



— m 2 + /i 2 ^ 2 , 

777/ + /Ujpr H : In 



ra l + fi 2 p 2 + Py'My _ A f l 
a/to 2 + /^p 2 



(59) 



Note that all the expressions ([571) - ([51?]) have finite limits at /i 2 — fi\ =0. From 
([5T)l - (|5"9"|) one can obtain: 

e - Pi, - 2P ± = 



(25 + ^toV 2 ^, 
2tt 2 



In 



m 2 + fijp 2 + pJ fif, - Mx \ 



a/to 2 + ll\p 2 



f( P 2 )pdp. (60) 



Besides, note that in the case fin ^ fi± it is essential that the logarithmic 
functions in the above formulae for the moments are to be substituted by the 
functions arcsin. 



4 Ultrarelativistic plasma 

4.1 Energy density and pressure 

magnetosphere (see, e.g., [H]), leads to very high values of the kinetic energy of 
electrons and protons (up to 10 18 ev). Thus GMSW appear to be always realized 
in an ultrarelativistic plasma. 
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In the ultrarelativistic plasma the particle rest mass does not affect the 
macroscopic moments in the drift approximation, therefore we can set: 



0||)e = 0||) P ; (p±)e = 
For an ultrarelativistic plasma, 
[i\_V ~> m ; n\\p > m ; 
the above expressions become explicit functions of and fi\\ : 



Pi = 



E fJ? ± »f\ 
Mi 



2 ni 



Mil 



MI 



My - Ml 



In 



/'II 



Mij - Mi 



M-L 



(61) 
(62) 

(63) 



P 



mImh 



4 Mm 



Mi 



2Mn ' Mi 



Mil 



Mi 



My - Mi 



M-L 



(64) 



^0 2 

e = -yMxMll 



where 



Mi 



Mij - Mi 



In 



OO 

^=E^ Jf(p 2 )P 3 dp = Const 



M|| + yMy - Mi 

M-L 



(65) 



(66) 



is the energy density of an isotropic homogeneous plasma. From (|60|) . in the 
approximation under study it follows: 



e - Pn - 2Pj = . 



(67) 



4.2 Special solutions 

With the expressions obtained, it is still very difficult to extract any information 
from Eqs. (|4"6"|) and (14T)1) . Therefore let us study these equations at the limiting 
values of the plasma anisotropy parameter: 



Mil 



(68) 



4.2.1 Isotropic plasma: 5 — >• 1 

Calcuating the limits of the expressions ((62 

P|| = Pi = U = pi 4 E„ , 

11 3 6 



66]) as <5 — > 1, we get: 



(69) 
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where (in = (i± = (i. Then it is easy to verify that Eqs. (pR))) and ([55]) have 
similar solutions: 

» 3 = §V¥- (70) 

This fact is of fundamental importance as it provides a microscopic explanation 
to RMHD pQ. In reality, if Eqs. gU) and ([56]) had independent consequences, 
it would indicate that the initial anisotropy of a collisionless plasma is broken 
by a gravitational wave, i.e. the RMHD equations would have been ineligible 
for the description of the collisionless plasma. 

Substituting relation ((70)) into Eq. we get the known solutions Q] for 

an isotropic ultrarelativistic plasma: 

* = ik:>> ^ 

e 2/3 

in this case 

M=£-4- (74) 



4.2.2 Plasma with a chilled transverse momentum: 6—^0 
In this case from (|62j) - (|66|) we obtain: 

1 

2 

Then Eq. (|46| has an integral: 



P„ « e = ^o/if /ii ; Px « . (75) 



e H . (76) 



Using in (|76| the integral (l56l) and the expression for H (jTTJ) , we obtain an 
expression for : 

M-L = (M-l)o— ■ (77) 

Now, substituting (1771) into (|75|) and the result obtained into (|46|) . using the 
integral (j5^|) . we get: 

e -2/3 

*' = Zw ; 

(A'||)o , 7a ^ 

" u = wT (78) 
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(79) 



From the solutions presented it is obvious that on the GMSW front (A — > 0) the 

o 

transverse pressure component is almost conserved (P± ~P±), whereas the the 
longitudinal pressure component tends to infinity proportionally to A^ 1 . thus, 
on the GMSW front the plasma initial anisotropy (in case it exists) intensifies. 

5 Anisotropy effect on the magneto - 
bremsstrahlung 

The squared projection values, averaged over the distribution ([29)) in the co- 
moving FR, are of the order 

< p\ >~ m± ; < p\ >~ /4 • (8i) 

therefore the average value of the ultrarelativistic particle energy is 

hl+pf\~ (82) 



thus in the case of a highly anisotropic plasma (for 6 — > 0) we obtain from ([82)) : 

< £ >~ M|| . (83) 
Substituting [i± and /xy from (l77|) and ([78)) into these expressions, we get: 

e -2/9 

< pi >«< pi > -p- « Const ; 

<£>«1^£. (84 ) 

lVa 

At small values of 5 the average value of the angle x between the vectors of 
microscopic velocity of an ultrarelativistic particle V and the magnetic field 
intensity H is small as well: 

SinX = T * ^ft^ =* * « Xoe^TA , (85) 
where xo = 5< 

As is known (see e.g. [12]), the magnetobremsstrahlung intensity for an ul- 
trarelativistic particle is concentrated in a narrow cone with the axis H and the 
angle 2% at the top. Thus, if 2\o is the output angle of magnetobremsstrahlung 
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of a nonperturbed plasma, then 2\ is the output angle of magnetobremsstrahlung 
in the GMSW. On the GMSW front (A(u) -> 0) the output angle of the radia- 
tion tends to zero. In this case the total intensity of the magnetobremsstrahlung 
of an electron jT2"]: 

according to ([501 and (jgijl in the GMSW is equal to: 

e -2/3(u) 

^-^l^WH (87) 

(where io is an unperturbed radiation intensity) and tends to infinity on the 
GMSW front. 

If the output angle of the magnetobremsstrahlung of the pulsar magneto- 
sphere is connected with the duration of the radiation pulse, then, as the grav- 
itational wave travels through the pulsar magnetosphere, a strong contraction 
of the pulse with a simultaneous increase of its intensity, will be observed. This 
effect offers an Alternative, as compared to [3] and [3J, opportunity to explain 
the giant pulses of the pulsars, whose radiation direction diagram is in most 
cases pencil-like. 
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